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Abstract 

The measurement of the deviation of the 2-3 leptonic mixing from maximal, D23 = 1/2 — sin^ 623, 
is one of the key issues for understanding the origin of the neutrino masses and mixing. In the Siy 
context we study the dependence of various observables in the atmospheric neutrinos on -D23- We 
perform the global 3i^-analysis of the atmospheric and reactor neutrino data taking into account 
the effects of both the oscillations driven by the "solar" parameters {Arn^i and 612) and the 1-3 
mixing. The departure from the one-dominant mass scale approximation results into the shift of 
the 2-3 mixing from maximal by Asin'^ 623 ~ 0.04, so that D23 ~ 0.04ib0.07 (Ifx). Though value of 
the shift is not statistically significant, the tendency is robust. The shift is induced by the excess of 
the e-like events in the sub-GeV sample. We show that future large scale water Cherenkov detectors 
can determine D23 with accuracy of a few percent, comparable with the sensitivity of future long 
baseline experiments. Moreover, the atmospheric neutrinos will provide unique information on the 
sign of the deviation (octant of 023)- 
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I. INTRODUCTION 



The present 2z/ analysis of the atmospheric neutrino results 0] terms of 
oscillations gives as mass squared difference and mixing: 



Am^2 = (1-3 - 3.0) X 10"^ eV^ sin^ 2^23 > 0.94 (90% C.L.) . (1) 

The results of SOUDAN and MACRO |5i experiments are in a good agreement with 
and this oscillation interpretation has been further confirmed by the K2K results 

The best fit of the data corresponds to maximal mixing sin^ 2^23 = 1-0, and this is 
one of the most striking results in neutrino physics. The maximal or close to maximal 
mixing implies a new symmetry of Nature which does not show up in other sectors of theory 
indicating its non-trivial realization. 

However, at the moment, one cannot claim that the mixing is indeed close to maximal. 
From the theoretical point of view the correct parameter which characterizes the deviation 
is not sin^ 26*23 but sin^ ^23 or 

D23 = ^-Sin2 023. (2) 

It is sin 6*23 that relates to the expansion parameter in the neutrino mass matrix and D23 
characterizes violation of the symmetry responsible for maximal mixing. From we obtain 
for the deviation 

ID23I < 0.12, (90% C.L.) , (3) 

and 

(O.5-sin2 023)/sin2^23 -0.3. (4) 

That is, the deviation can be of the order of mixing itself. 

In the case of a significant deviation we cannot speak of a special symmetry, and in fact, 
the large 2-3 mixing may appear as the sum of small (order of Cabibbo) mixing angles (see, 
e.g., So, maximal or non-maximal mixing is equivalent to the dilemma of new symmetry 
or no-new symmetry of Nature. (Here we exclude the possibility that small angles sum up to 
give accidentally exact maximal mixing). Depending on value of the deviation the approach 
to the underlying physics can be different. 

The present data may already give some hint of deviation of the 2-3 mixing from maximal. 
Indeed, there is some excess of the e— like events in the sub-GeV range. The excess increases 
with decrease of energy within the sample 0| . In comparison with predictions based on the 
atmospheric neutrino flux from Honda |^ the excess is about (12-15)% in the low energy 
part of the sub-GeV sample (p < 0.4 GeV, where p is the momentum of lepton) and it 
shows no significant zenith angle dependence. In the higher energy part of the sub-GeV 
sample (p > 0.4 GeV) the excess is about 5%, and there is no excess in the multi-GeV 
region {p > 1.33 GeV). 

In principle, the observed excess is within the estimated 20% uncertainty of the original 
atmospheric neutrino flux. So the 2z/ analysis of data with free overall normalization and 
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tilt of the energy spectrum can explain a large enough fraction of the excess leading to the 
result of best fit maximal 2-3 mixing. 

The excess has become more significant in the latest Super-Kamiokande analysis The 



recent data on primary cosmic rays 
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as well as the new 3-dimensional calculations 



of the atmospheric neutrino fluxes 1^ imply a lower neutrino flux, and therefore a larger 
excess which is becoming more difficult to explain by a change of the atmospheric neutrino 
fluxes within their expected uncertainties 0, |^ . 

Alternatively, such an excess can be explained (at least partly) by the i/e-oscillations 
driven by the solar oscillation parameters [l3[ provided that the 2-3 mixing deviates from 
maximal Q|. For the solar parameters which we will call the LMA parameters the 
combined analysis of the solar and KamLAND [l^ data leads to 3a ranges [f^ : 



Am2i 



(7.4 - 9.2) X 10"^ eV^ 



tan^ 9i2 = 0.28 - 0.58 , 



(5) 



Oscillations of atmospheric neutrinos drive n by these LMA parameters have been widely 
discussed in the literature 
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23| . It is found that the relative change 



of the atmospheric z/g flux due to oscillations driven by the solar parameters is determined 
by the two neutrino transition probability P2(Am2i, 6*12) and a "screening" factor 1^ : 



- 1 = P2{Aml„9,2) (rcos2 023 - 1) 



(6) 



where and are the electron neutrino fluxes with and without oscillations, and r = 
F^/F^ is the ratio of the original muon and electron neutrino fluxes. The screening factor 
(in brackets) is related to the existence of both electron and muon neutrinos in the original 
atmospheric neutrino flux. 

In the sub-GeV region r 2, so that the screening factor is very small when the — Vr 
mixing is maximal. According to Eq. (jHI), the excess of the e-like events can be written as: 



(7) 



where {P2)uv = [(1 ~ Oi^^) + ^(A)], and (P2) ((A)) is the average transition probability 
for neutrinos (antineutrinos) in the Earth matter. The parameter ^ gives the relative contri- 
bution of antineutrinos (without oscillations). For the sub-GeV electrons we have ^ ~ 0.3. 
Once the solar oscillation parameters have been well determined, one can calculate Po ra ther 
precisely. Then the study of the excess can be used to search for the deviation D23 |lq . 

For the presently allowed range of solar oscillation parameters, neutrino oscillations can 
lead up to a (5-6)% excess of the e-like events in the sub-GeV atmospheric neutrino sam- 
ple Q Q. So, the oscillation explanation of the observed excess would imply maximal 
allowed deviation 1)23 • On the other hand, the decrease of sin^ 2^23 influences other observ- 
ables (like high statistics measurements of the zenith angle dependence of the //-like events). 
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Therefore to make a definitive conclusion about the deviation one needs to perform a com- 
bined analysis of the all available data and to take carefully into account the uncertainties 
in the atmospheric neutrino fluxes. 

The z/g oscillations are also induced by non-zero 1-3 mixing and Amg^ responsible for 
the dominant mode of the atmospheric neutrino oscillations. This effect is mostly visible 
for the multi-GeV sample 3 25. bd] for which the Earth matter effect becomes important 
and can enhance the oscillations. Non-zero 1-3 mixing induces also an interference effect in 

iowever, within the present bound on the 1-3 mixing from the 
irt - the dominant factor which leads to a possible excess of the 



the sub-GeV range [H 0- 
CHOOZ reactor experiment 



sub-GeV e-like events is the Am^^-driven transitions discussed here. 

In this paper we perform a detailed study of the dependence of the atmospheric neutrino 
observables on the deviation D23. We determine D23 from the analysis of present data and 
investigate possibilities of future experiments. 

The paper is organized as follows. In Sec. [TTl we discuss the dependence of different 
samples of the atmospheric neutrino data on D23- In Sec. II I II we present the results of the 
global analysis of the atmospheric and CHOOZ results in terms of three-neutrino oscillations 
where the effect of both mass differences is explicitly considered. In Sec. IIVI we study the 
capabilities of future large scale water Cherenkov detectors to determine D2z- Discussion of 
the results and conclusions are given in Sec. El 



II. i:>23 AND THE ATMOSPHERIC NEUTRINO OBSERVABLES 

In this section, using some approximate analytical results, we discuss the dependence of 
the atmospheric neutrino observables on the deviation D23- 

1. According to Eq. (|7j) the excess of e-like events in the sub-GeV range due to LMA 
parameter oscillations is proportional to the deviation D23- 

ee ^ D2; r (P2).p, (8) 

while their zenith angle distribution (encoded in (^2)^^) does not depend on D23- 
The 1-3 mixing modifies this dependence but, for values compatible with the CHOOZ 
bound, it is a subdominant effect for sub-GeV events. 

The important point is that the excess decreases with energy as ~ E'"^. This par- 
ticular energy dependence of the excess allows to disentangle it from the uncertainties 
of the neutrino fluxes. 

2. As a consequence of this energy dependence, in the multi-GeV range the excess of the 
e-like events due to the LMA parameters is very small: 5-10 times smaller than in 
the sub-GeV range and therefore below 1%. Conversely, the zenith angle distribution 
here is stronger. 
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3. The sub-GeV //-like events have more comphcated dependence on the deviation [3 



N, 



- 1 = - sin^ 26'23(sin^ -)uu-e^- ■ (9) 

The dominant contribution due to the z/^ — Uj. oscillations, P^,-, (the first term on the 
RHS) depends on sin^ 26*23 = 1 ~ ^Dga- Here is the oscillation phase due to Aml^. 
The term describes the decrease of the rate of the /x-like events due to oscillations 
of the muon neutrinos into the electron neutrinos driven by the 1-2 mixing: 

~ ^^^ee ^ 1^23 COS^ ^23(P2).P • (10) 

The third term is the interference of these two contributions and it is essentially 
averaged out. The expression can be rewritten as 

^ - 1 ^ {Wl, - l)(sin2 - D23 COS^ 923{P2)u, . (11) 

Notice that the 1-2 mixing has an opposite effect on the rate of the /x-like events as 
compared to its effect on the e-like events. Thus for the case of the excess of e-like 
events, the 1-2 mixing diminishes the rate of /i-like events. So, it cancels partially the 
increase of the rate due to the main term P^^. Furthermore, both terms exhibit a 
different dependence on D23: 

P^r^il-^Dl^), eM(xP)23. (12) 

So, for large deviation the change of the main term dominates, whereas for small 
deviations the two contributions become comparable. 

4. For muons in the multi-GeV range, due to the suppression of P2, the effect of the 
1-2 mixing is small and the change of the main term in Eq. Q dominates down 
to very small deviations. So, basically the rate of the /z-like events increases with 
|D23|- Conversely, for the multi-GeV events the effect of 1-3 mixing can be more 



substantial 



mm, 



5. In the sub-GeV range the double ratio can be written as Q| 



n^/e - n^/e 7— , (i-sj 

-L T te 

where P™^^ is the double ratio in the case of two-neutrino oscillations with maximal 
mixing. In terms of the deviation it can be rewritten as 

o pmax 0-5 + 2Z^i3-^23COS^g23(P2).P .... 

"^'^^ - "^-/^ l + D23r{P2U ' ^ ^ 
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For D23 > the 1-2 mixing effects partially compensate the change of main term and 
numerical inspection of Eq. (jl4|) reveals that the change of -R^/e with D23 is rather 
weak. For D23 < the 1-2 mixing enhances the ratio R^/e- 

In the multi-GeV range the e corrections are small and the double ratio increases with 
D23. Comparing the double ratios in the sub-GeV and multi-GeV ranges we conclude 
that -R^/e(sub — G) changes weakly (for D23 > 0) whereas i?^/e(multi — G) increases, 
so that the ratio of ratios 

/?^/e(multi-G)/i?^/e(sub-G) (15) 

increases with |-D23|- 

6. For upward-going muons the average energy of the neutrinos is above 10 GeV. For these 
energies, the effect of 1-2 mixing is strongly suppressed in matter. Also the possible 
1-3 mixing has additional matter suppression: the relevant factor is 2EV / /S.im? ~ 0.1. 
So, to a good approximation the rate depends on 1 — 4Z)23- 

Summarizing, the rate of the low energy e-like events is proportional to the deviation 
D23 and the rates of high energy (/x-like) events depend on (1 — 4D23)) ^-fi-; the deviation 
squared. The rate of low energy /x-like events an the double ratios may have non-trivial 
interplay of the two dependences: cancellation or enhancement of the main mode and 1- 
2 mixing contributions. With these consideration in mind one can better understand the 
results of the various analyses described in the following sections. 



III. GLOBAL ANALYSIS OF PRESENT DATA 

We discuss here what can be learned from the subleading effects induced by non-vanishing 
solar splitting Am|^ and the present atmospheric and reactor neutrino data. In our analysis 
we include the complete 1489-day charged-current data set for SK-I which comprises the 
sub-GeV and multi-GeV e-like and //-like contained event samples (each grouped into 10 
bins in zenith angle), as well as the stopping (5 angular bins) and through-going (10 angular 
bins) up-going muon data events. In the calculation of the event rates we have used the new 
three-dimensional atmospheric neutrino fluxes given in Ref. jl^ . 

Details of our statistical analysis based on the pull method can be found in the Appendix 
of Ref. and here we summarize some points which are essential for the present study. 
Together with the statistical errors, we consider two types of uncertainties: the theoretical 
and systematic ones. 

The theoretical uncertainties include uncertainties in the original atmospheric neutrino 
fluxes and in the cross-sections. 

We have parametrized uncertainties of the atmospheric neutrino fluxes in terms of four 
pulls: 
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FIG. 1: Allowed regions (90%, 95%, 99% and Scj C.L.) of the oscillation parameters Am^g and 
sin^ 023 from the analysis of different atmospheric data samples. The best fit points are marked 
with either a star (Arrigi ^ 0) or a triangle (Am^x = 0) . Colored regions and stars correspond to 
Am^i = 8.2 X 10~^ eV^ and tan^ 9i2 = 0.42, whereas hollow regions and triangles are for Am2i = 0. 
In both cases we assume ^13 = 0. 



a total normalization error, which we set to 20%; 

a "tilt" factor which parametrizes possible deviations of the energy dependence of the 
atmospheric fluxes from the simple power law defined as 



E 



ME) 



(5 In 



E 

En 



(16) 



with an uncertainty on the factor 6, erg = 5% and Eq = 2 GeV; 

the uncertainty on the i/^/ z/e ratio, which is assumed to be a^/e = 5%; and 
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• the uncertainty on the zenith angle dependence which induces an error in the up/down 
asymmetry of events which we conservatively take to be 5%. 

We also include independent normalization errors for the different contributions to 
the interaction cross section: quasi-elastic scattering (QE), cxnorm = 15%, single pion 
production(l7r), aj^orm = 15%, and deep inelastic (DIS) scattering (also refer to as multi-pion) 
for which we estimate cXnorm = 15% for contained events and cr'^oTm = upward-going 
muons.^ 

We include as systematic uncertainties the experimental uncertainties associated with 
the simulation of the hadronic interactions, the particle identification procedure, the ring- 
counting procedure, the fiducial volume determination, the energy calibration, the relative 
normalization between partially-contained and fully-contained events, the track reconstruc- 
tion of upgoing muons, the detection efficiency of upgoing muons, and the stopping-thrugoing 
separation. 

In order to illustrate which data samples are more sensitive to the departure from the one- 
mass-scale dominance approximation and to the deviation D23 we first perform the analysis 
for different sub-samples. The results of these partial analysis are presented in Fig. ^ where 
we show the allowed regions in the {Amli, sin^6'23) plane. The colored (shadowed) regions 
correspond to Am|^ = 8.2 x 10~^ eV^ and tan^ 612 = 0.42, whereas hollow regions are 
for Amli = 0. In both cases we assume ^13 = 0. A comparison between the two sets of 
regions clearly shows that the main effect of Amgi oscillations appears in the e-like events 
at lower energies as discussed in the previous section. As seen in the figure the inclusion of 
Am2i-driven oscillations in the analysis breaks the symmetry of 6*23 around maximal mixing 
providing the expected sensitivity to D23- In accordance with the considerations of Sec. |n] 
the allowed regions and the best fit point shift to sin^ 6*23 < 1/2. 

Also, as discussed in the previous section, the effect is much more suppressed at higher 
energies. As can be seen in the figure, the high energy muon neutrino events as well as the 
multi-GeV (e-like and /i-like) events do not lead to the shift of 2-3 mixing from maximal. The 
contained /i-like events produce a very small shift in the opposite direction: to sin^ 6'23 > 1/2. 
Let us underline that despite they show much less sensitivity to 1-2 oscillations, the muon 
neutrino data are very important to constraint the 2-3 oscillation parameters and to limit 
the effect of theoretical and systematic uncertainties. 

Next we study the effect of Arn^i oscillation in the combined analysis of all available 
atmospheric neutrino data (a total of 55 data points). In order to account also for the 
effect of the angle ^13 we include in the analysis the results of the CHOOZ experiment. 



For CHOOZ we consider the energy binned data. T, 
positron spectra from both reactors. Table 4 in Ref. 



lis corresponds to 14 data points (7-bin 
) with one constrained normalization 



^ We also account for the uncertainty of the ai^v^J(^i,v^ ratio which is relevant only for contained events, 
and it is much smaller than the total normalization uncertainty. 
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FIG. 2: Allowed regions (90%, 95%, 99% and 3a C.L.) of the oscillation parameters Am|;^ and 
sin^ 023 from the combined analysis of all the atmospheric and CHOOZ data samples. The best fit 
points are marked with either a star or a triangle. In the lower panels we show the dependence of 
the function on 623, marginalized with respect to Am|^. Colored regions, stars and solid blue 
lines correspond to tan^ 612 = 0.42, ^13 free and Am^]^ set to the value indicated in each panel. 
Hollow regions, triangles and dashed black lines are for Am^i = and 6*13 free. 



parameter. In this analysis we have assumed CP conservation but we have considered both 
possible values of the CP parity which correspond to the CP phases 6 = and 6 = n. 

The results are shown in Fig. |21 As before, the colored (shadowed) regions correspond to 
tan^ 612 = 0.42 and Amgi different from zero, whereas the hollow regions are for Arrigi = 0. 
In order to verify explicitly that our results are robust and do not change for non-zero ^13, 
we have marginalized Xsk+chooz with respect to this parameter. In the lower panels we 
plot the function marginalized with respect to Arrig^^ as well. 

From the figure we see that, even with the present uncertainties, the atmospheric data 
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has some sensitivity to Am^i oscillation effects and that these effects break the symmetry 
in 623 around maximal mixing, although the effect is small. Quantitatively, for Amgi = 
8.2 X 10~^ eV^ and tan^ 6^12 = 0.42, the best-fit point is located at sin^ 6^23 = 0.46, with la 
(3a) interval (0.30) 0.39 < 8111^623 < 0.54 (0.65) which means that for these values of the 
solar parameters we find 



whereas for Amli = we obtain D23 = 0.0 ± 0.07 (^017) (the slight asymmetry of errors 
in this case is induced by ^^13). For values of Amgi in the range indicated by solar and 
KamLAND data, the deviation from maximal mixing of ^23 increases with Am2i. We have 
also verified that, once CHOOZ is also included in the analysis, it makes little difference to 
leave ^13 free or to set it to zero. 

Let us stress that, although statistically not very significant, this preference for non- 
maximal 2-3 mixing is a physical effect on the present neutrino data, induced by the fact 
than an excess of events is observed in sub-GeV electrons but not in sub-GeV muons nor, 
in the same amount, in the multi-GeV electrons. As a consequence, this excess cannot be 
fully explained by a combination of a global rescaling and a "tilt", of the fluxes within 
the assumed uncertainties. In the pull approach we find that both the total normalization 
and the tilt pulls are essentially fixed by the combination of low-energy and high-energy 
muon data, and there is no freedom left to accommodate the remaining excess of low energy 
electron events. Such an excess can only be partially explained by means of another pull, the 
/i/e flavor ratio, whose uncertainty, 5%, is however much smaller. We have explicitly verified 
that in the vicinity of the best-fit point the only pull which is affected by the precise value 
of D23 is the fi/e flavor ratio, whereas the total normalization and the tilt are practically 
insensitive to it. In particular, we have checked that increasing the tilt uncertainty by a 
factor of 2 or allowing for a totally unconstrained overall normalization does not affect the 
present value of D23 in Eq. (fTTjl . 

From the lower panels in Fig. |21we can also see that the quality of the fit slightly improves 
when Am^i differs from zero. As expected, this is due to the fact that the non-vanishing 
value of Am2i and non-zero D23 imply that z/^ v^. transition is more efficient than the 
inverse one, Vf, z/^, which allows to partially explain the excess of e-like events observed 
by Super-Kamiokande in the sub-GeV data sample. 

Notice that the central value of deviation in Eq. (fTTj) corresponds to sin^ 26^23 = 0.9936 
which is beyond the sensitivity of the next generation of the long-baseline experiments. 

In summary, in this section we have shown that atmospheric neutrino data are sensitive 
to the subleading z/g conversion induced by a non vanishing IS.m^^. More important, 

the preference of atmospheric data for maximal ^23 mixing appears to be a specific property 
of the one-dominant mass scale approximation, and seems to disappear when oscillations 
with the two wavelengths between all three known neutrino flavor are considered. However, 
present data still have far too little statistics to provide a conclusive answer. 




(17) 
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IV. SENSITIVITY OF FUTURE EXPERIMENTS 



Having shown that atmospheric neutrino data can be a useful instrument to search for 
deviations of ^23 from 45°, we now discuss what can be learned from future atmospheric 
experiments. For the sake of concreteness we have assumed a SK-like detector with either 
20 (SKx20) or 50 (SKx50) times the present SK-I statistics and the same systematics as 
SK-I, and we have used the same event samples as in SK. 

The procedure is as follows: First we simulate the signal according to the expectations 
from some specific choice of the "true" values of parameters which we denote by aJ 



and then we construct 



uj = (Am^i, Am^^, ^12, ^13, ^23) , (18) 



XL(A"^21> A"^31> ^12, ^13, ^23 I ^) (19) 



assuming 20 or 50 times the present SK statistics and three choices for the theoretical and 
systematic errors (see definitions in Sec. Illljl : 

(A) same theoretical and systematic errors as in present SK; 

(B) same systematic errors as in present SK, but no theoretical uncertainties; 

(C) neither theoretical nor systematic uncertainties (perfect experiment). 

Next, in order to study the effect that non-zero values of Am^i and 612 can produce in the 
determination of the atmospheric parameters Am|^ and ^23 we define 



xL™(Am^„e23|^)- min 



Xsk(A"^21' A"^31' ^12 = ^12, On, \ Lu) 

+ XcHooz(A"^2i, ^rnli. ^12 = O12, Oi3 \uj)+' ^""21 - A'-2i 



2- 

(20) 



where we minimize with respect to the solar and reactor parameters Am^^ and 613 and we 



^21 
,2 



keep only the explicit dependence on the "atmospheric" parameters Am^i and 623. The as- 
sumption 612 = 612 is made for purely practical reasons because a complete scan of the whole 
five-dimensional parameter space requires too much computer time. Note that regardless of 
the specific assumptions on the 'true values' Am^i and 6^13 the parameters Am^i and 6*13 
are allowed to vary in our fit. In the definition of Xatm+reac ^'i- we have included 
also the CHOOZ experiment Xchooz order to have a realistic bound on ^13. Similarly, the 
term [(Am2x — Am^i) / '^Aml^)]'^ accounts for the bound on Am^i which is expected from 
KamLAND in the next few years. Following Ref. we have assumed that by then Am2i 
will be known with an uncertainty of 3% at la. 
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FIG. 3: Allowed regions (at 90%, 95%, 99% and 3a C.L.) of oscillation parameters Am^g and 
sin^ 023 expected from an atmospheric neutrino experiment with 20 times the present SK statistics 
and the same theoretical and systematic errors as in present SK. For definiteness, we choose 
013 = 0, Am2i = 8.2 X 10~^ eV^ and tan^ 0i2 = 0.42, and we scan different values of Am|^ and 023- 
We also include the constraints from the CHOOZ experiment, as well as the sensitivity to Am^i 
expected after 3 years of KamLAND data (Eq. (|20|) ). The undisplayed parameters Am^x and ^13 
are marginalized. 
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As an illustration of the expected sensitivity from future atmospheric neutrino experi- 
ments, we show in Fig.|21the allowed regions obtained from Xatm+reac assuming 20 times the 
present SK statistics and the same theoretical and systematic errors as in present SK (case 
A). For definiteness, we choose 6'i3 = 0, Am^i = 8.2 x 10~^ eV^ and tan^ 9i2 = 0.42, and we 
scan different values of Amg^ and ^23- From this figure we find that: 

• future atmospheric neutrino experiments can observe and measure deviations of 623 
from maximal mixing, provided that 623 is not too close to 45°: sin^ ^23 < 0.38 or 
sin^ 6*23 > 0.60; future reduction in the theoretical errors will further improve the 
sensitivity; 

• they can discriminate between the "light-side" and "dark-side" for 6*23, i.e., they are 
sensitive to the octant of 6'23. 

In the rest of this section we quantify these two possibilities. 



A. Deviations from maximal mixing 



Let us quantify the sensitivity of the future atmospheric neutrino experiments to deviation 
of ^23 from 45°. We compare our results with the corresponding bounds which can be 



expected for future long baseline (LBL) experiments 29|, |30|, following the discussion in 
Ref. [29]. The results of our analysis are summarized in Fig. 0]and Table HI which can be 
directly compared with Fig. 1 and Table 1 of Ref. ji^. 



To perform this analysis, we have constructed the following function: 



Ay' 



mm 



ATM + RE AC 



(Am^i, ^23 = 45°|oJ)-x 



ATM+REAC 



(Amgi, 6123 \uj) 



(21) 



ATM+REAC (^"^31' ^23 1"^) is given in Eq. ((201). Fig. IH we plot the dependence of 



where xl 

^Xno-max ^"^31 ^23, for both Am2i = 8.2 X 10"^ eV' (colored regions) and Am2i = 
(hollow regions). We take tan^ 6'i2 = 0.42 and ^13 = 0. The blue, green and yellow regions 
correspond to Axno-max > I5 4 and 9, respectively. In other words, in Fig. HI we display, for 
each value of Aml^, the range of 6'23 for which the simulated signal can be reconstructed 
as having maximal ^23 at 1, 2 and 3 a. The white region corresponds to the the range of 
6'23 for which the simulated signal cannot be distinguished from maximal 6'23 at Icr. The 
corresponding bounds on 6'23 for Am|^ = 2.2 x 10~^ eV' are summarized in Table HI 

From Fig. |3] and Table IJ we find that the sensitivity of atmospheric neutrino data to 
deviations from maximal mixing for large values of Amg^^ is comparable to what can be 
expected "after ten years" from LBL experiments according to Ref. 29|, D23 < 0.050 (0.069) 



at 90% (3cr) CL. Furthermore, for small values of Am|]^ the atmospheric neutrino studies are 
much more sensitive than LBL experiments, which lose sensitivity very fast when Amg^ < 
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FIG. 4: (Am|]^, ^23) regions with Axno-max smaller than 1 (white), between 1 and 4 (blue), between 
4 and 9 (green), and larger than 9 (blue), respectively for Am2^ = 8.2x 10~^ eV^ (the hollow regions 
are obtained for Am|^ = 0). We set ^13 = and tan^ ^12 = 0.42. The dashed oval marks the 3a 
region presently preferred by SK data. 



Experiment 






[0.5 - sin2 


^23 








With Am^i 


= 8.2 X 


10^5 eV^ 


With Am|i = 






90% C.L. 


3cr 


90% C.L. 


3cr 


SKx20 (A) Theo+Sys+Stat 


[-0.086, 0.067] 


[-0.116, 


0.096] 


-0.080, 0.080] 


[-0.108, 


0.108] 


SKx20 (B) Sys+Stat 


[-0.040, 0.050] 


[-0.062, 


0.075] 


-0.058, 0.058] 


[-0.078, 


0.078] 


SKx20 (C) Stat only 


[-0.032, 0.032] 


[-0.054, 


0.052] 


-0.054, 0.054] 


[-0.073, 


0.073] 


SKx50 (A) Theo+Sys+Stat 


[-0.070, 0.053] 


[-0.094, 


0.077] 


-0.064, 0.064] 


[-0.087, 


0.087] 


SKx50 (B) Sys+Stat 


[-0.030, 0.040] 


[-0.046, 


0.061] 


-0.046, 0.046] 


[-0.063, 


0.063] 


SKx50 (C) Stat only 


[-0.021, 0.021] 


[-0.037, 


0.036] 


-0.042, 0.042] 


[-0.058, 


0.058] 



TABLE L Rejection of maximal mixing expected from future atmospheric neutrino experiments. 
We assume ^13 = 0, tan^ Ou = 0.42 and Am|^ = 2.2 x 10~^ eV^, and we study both the case 
Am^i = 8.2 X 10"^ eV^ and Am^^ = 0. The given intervals correspond to Axno-max (see Eq. I\'21\i ) 
smaller than 2.71 (90% C.L.) and 9 {3a). 
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2 X 10~^ eV^ while the bound which can be obtained from the atmospheric neutrino data is 
practically independent of the value of Am|j^ . 

The comparison among the left, central and right panels of Fig. |3] also shows that the 
sensitivity of atmospheric neutrino data to deviations from maximal mixing improves con- 
siderably if theoretical errors on the atmospheric fluxes and cross sections are reduced. On 
the contrary setting to zero the systematic uncertainties induce a smaller improvement. This 
implies that the obtained results hold even if the future atmospheric neutrino experiment is 
affected by somewhat larger systematics than the present SK detector has. 

We also see that, as expected, when Arnl^ ^ the ranges of 623 can be asymmetric. This 
effect is mostly seen in the first two panels (cases A and B) because larger errors allow for 
larger values of D23. We find that the overall effect of the theoretical errors is such that 
the fit for maximal mixing is "less-bad" if an excess of e-like sub-GeV events is observed 
as compared to the observation of a deficit, while for systematic uncertainties the opposite 
holds. 

In any case, comparing the solid (obtained with Afn^i = 8.2 x 10~^ eV^) and the hollow 
(obtained with Am^i = 0) regions in Fig. I^we see that the value of the solar mass splitting 
is not the most important effect in the discrimination from maximal mixing, and the bound 
comes mainly from muon data. Only when both theoretical and systematic uncertainties 
are neglected (case C) the bound on D23 becomes visibly sensitive to Am^i- This occurs 
because the effect of a non-zero value of Am^i is comparable to the small statistical error so 
this small effect is relevant only when the fit is purely statistics-dominated. In general, the 
subdominant Am^i effect is mostly important to determine the octant of 623 as we discuss 
next. 

B. Determination of the octant 

As an illustration of the capability of atmospheric neutrino data to discriminate between 
023 smaller or larger than 45°, we show in Fig.^lthe dependence of Xatm+reac ^ particular 
simulated point U as a function of 623, after marginalizing over all other parameters. In what 
follows we work under the hypothesis that by the time this future atmospheric neutrino 
experiment is in place we have found no-evidence of 613 but we may have a better-than- 
present determination of the oscillation parameters from terrestrial experiments. To account 
for this effect we have constructed the function 




,2 



2 



,2 



2 



2 



+ XcHOOz(^"^2l5 ^"^31' ^12 = O12, 9l3 I ^) + 




15 




FIG. 5: Dependence of Xatm+reac (^^1^ black line) and Xatm+reac+lbl ('^^'^ dashed line) on sin^ 623, 
for Aml-^ = 2.2 x 10"'^ eV^ and setting the simulated point (uJ) to 6*13 = 0, tan^ 612 = 0.42, 
Am|i = 8.2 X 10"^ eV^ sin^^as = 0.42 and Am^^ = 2.2 x 10"^ eV^ For the definition of 
xixM+REAc and xiTM+REAc+LBL scc Eqs. (jSni) and ((22I). 



+ 



/ Amj-^ - Am^^ V + ^^^^ ^^23 - sin^ 2^23 \ ^( ^^^^ 2^13 - sin^ 2^13 ^ ^ 

\ ^^ml^ j \ ^sin2 2e23 / \ sin' 29^3 



(22) 



where in addition to the CHOOZ and KamLAND-3yr bounds we have included a stronger 
bound on the Oi^ angle (for example, from some future reactor experiment js^) as well 
as an improved measurement of the atmospheric parameters Am|]^ and ^23 from future 
narrow beam LBL experiments such as T2K [3^ or NuMi js^l- Following Ref. [s^, we have 
assumed crgjn2 251^3 = 0.01, (^^^^220^^ = 0.015 and a^jn2^/ Afnl^ = 0.015. Note that in Eq. 
we have expressed the sensitivity of future LBL experiments in terms of sin^2^23; rather 
than sin^6'23, to account for the fact that these experiments have no sensitivity to the ^23 
octant if ^13 turns out to be very small {i.e., they cannot distinguish between sin^ ^23 < 0.5 
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FIG. 6: Dependence of Ax^c (^ee Eq. on 623, for 6*13 = 0, tan^ 6I12 = 0.42, Am|i = 

8.2 X 10~^ eV^ and Amg^ = 2.2 x 10^^ eV^. The different lines correspond to the same cases as in 
Fig.ini 



and sin^ 623 > 0.5)^ j29 



In Fig. El we see that the best fit point (x^ = by construction) is located at sin^ ^23"'^ = 
sin^^23 (and Amg^ = Afnli, Arrigi = Am2i, and 613 = 613), while in the left panels (case 
A) presents a second local minimum at sin^ ^23'^^ ~ 1 — sin^^23- The shift between the 
position of the secondary minimum and the mirror symmetric value of the true minimum 
is expected from the atmospheric analysis. On the other hand, in the central (case B) and 
right (case C) panels the second "false" minimum has disappeared. Again, this illustrates the 
importance of improving our knowledge of the fluxes and cross-sections. As before, for both 
SKx20 and SKx50 the dominant source of errors is the theoretical uncertainties, whereas 
systematic uncertainties play a somewhat smaller role. We find that once the theoretical 
uncertainties are neglected, the atmospheric data can totally lift the degeneracy between 



^ A very long baseline wide beam experiment such as the BNL proposal could be also sensitive to 



subdominant Am^i effects. 
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the two octants of ^23- 

In order to quantify the discrimination power of the octant of ^23 we define the difference 









— min 




Am2j L 



^Xdisc(^) = min XATM+REAC(+LBL)(^"^31.^2f'' "mill XATM+REAC(+LBL)(^"^31.^2r 



UJ 



(23) 

where 6^^^"^ is either the mixing angle of the secondary local minimum, or (90° — ^23"*^) 
there is not a secondary local minimum. In Fig. El we plot this difference as a function 
of the simulated "true" value sin^^23 for tan2^i2 = 0.42, Am^^ = 8.2 x 10"^ eV^ and 
Aml^ = 2.2 X 10^3 eV^ This function is very insensitive to the exact value of Am|| in the 
interval 10~3 eV^ < Am|i < 5 x lO^^ eVl 

The figure shows that unless 6'23 is very close to maximal mixing such a future atmospheric 
neutrino experiment can provide a meaningful determination of the octant of 623. For 
example, from the figure we read that for SKx50 the octant of ^23 can be determined at 
90% CL if 

sin^ ^23 < 0.42 [023 < 40°] or sin^ ^23 > 0.57 [^23 > 49°] (A) 

sin^ 023 < 0.48 [^23 < 43°] or sin^ 623 > 0.52 [^23 < 46°] (B) (24) 

sin^ ^23 < 0.49 [^23 < 44.4°] or sin^ ^23 > 0.51 [^23 < 45.6°] (C) 

These results are almost independent of the exact value of Amg^ within the present atmo- 
spheric region. From Fig. IHl we see that the discriminating power can be slightly improved 
if LBL experiments provide a better determination of I-D23I, as a consequence of the shift in 
the position of the secondary minimum in the atmospheric neutrino analysis. 

A final comment on the role of ^13. If ^13 is not very small, the oscillation probabilities at 
the considered LBL experiments are also not symmetric under the change of octant and they 
can also contribute to the octant discrimination j^^. This effect has not been statistically 
quantified in detail in the literature and it is beyond the purpose of this paper. 



V. CONCLUSION 

In this paper we have discussed the phenomenology of atmospheric neutrinos associated 
with the deviation of the 2-3 leptonic mixing from maximal. Our main results can be 
summarized as follows: 

1. We have performed the global 3z/-analysis of the atmospheric and reactor neutrino data 
taking into account the effect of both the oscillations driven by the "solar" parameters 
{Am2i and 6*12) and of the 1-3 mixing. The results are shown in Fig. |21 We find 
that the departure from the one-dominant mass scale approximation in the analysis 
results into the shift of the 2-3 mixing from maximal, so that D23 ~ 0.04 ± 0.07. 
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The shift is due to the excess of e-hke events in the sub-GeV sample as illustrated in 
Fig. n For these values of D23 the LMA oscillations explain the excess only partly. 
Larger deviation is disfavored by the zenith angle distribution of the /i-like events. 
The qualitative effect of the shift of 2-3 mixing from maximal one is robust. Though 
particular value of the shift depends on details of the treatment of errors. 

2. Future experiments will have much higher sensitivity to D2z- With 20 — 50 SK statis- 
tics and better knowledge of the cross-sections and the original fluxes the atmospheric 
neutrinos will probe D23 down to few percent (see Fig. ^ and Table E} - a sensitivity 
comparable with that attainable at future LBL experiments. This sensitivity does not 
change with decrease of Am|^ and therefore a high statistics atmospheric neutrino ex- 
periment is better than LBL experiments to determine deviations of 6*23 from maximal 
mixing if Amg^^ lies in the lower part of the present allowed range. 

3. If 1^23 7^ future atmospheric neutrino experiments have the potentiality to discrimi- 
nate the octant due to effects associated to the LMA oscillations, as shown in Figs. 
and El 
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